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ff^ I Abstract 

We use exact WKB analysis to derive some concrete formulae in singular quan- 
tum perturbation theory, for Schrodinger eigenvalue problems on the real line with 
polynomial potentials of the form (g*^ + gq^), where > M > even, and 5 > 0. 
Mainly, we establish the limiting forms of global spectral functions such as 

the zeta-regularized determinants and some spectral zeta functions. 



1 Introduction 

The purpose of this work is to set up a path to obtain precise statements of a quantum 
perturbative nature with the help of exact WKB analysis. The RIMS has always played a 
major and pioneering role in the inception and growth of exact asymptotic analysis, and 
earlier, in the development of some of its fundamental tools (such as hyp erf unctions and 
holomorphic microlocal analysis). This influence is testified by the Proceedings volume 
of a recent Kyoto conference, which offers a very complete view of the subject 0. It is 
therefore a great honor and a proper tribute to RIMS to write about exact WKB analysis 
in this anniversary issue. 

A prototype problem in quantum perturbation theory is the quartic anharmonic os- 
cillator, 

/ d^ \ 

[-j^ + q^ + gq^-E)^>{q) = Q, g G M, (7 > 0. (1) 

This problem has a purely discrete eigenvalue spectrum {Ek{g)} for all (7 > 0. A typical 
task in (Rayleigh-Schrodinger) perturbation theory is to compute individual eigenvalues 
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Ek{g) (or their eigenfunctions) as formal power series of the couphng constant g pp. This 
is practically important when the unperturbed {g = 0) problem is exactly solvable, here a 
harmonic oscillator; a major drawback is however that the coupling term has the higher 
degree, hence the formalism is singular. Thus, the perturbation series converges for no 
(? 7^ 0; it only gives an asymptotic expansion for g —>■ 0, moreover non- uniform in the 
quantum number k. 

As our theoretical discussion can readily include all binomials potentials, we will ac- 
tually study the more general Schrodinger equation 

(-^ + f/p(g)-^)^(g) = 0, f/g(g) = + g e M, iV > M > 2 even, g>0; 

(2) 

we keep Ek{g) as a generic notation for the eigenvalues of this problem [{N, M)-dependences 
now being implied). 

(Exact WKB formalisms accommodate non-even potentials as well |21[Z1; for instance, 
eq.Q could be considered with odd or M but on the half-line [0, -|-oo) [Zj; however, 
this extension is not essential here while it does complicate the classification when M = 1, 
so we omit it in the present work.) 

A very basic fact {Symanzik scaling property) is that a simple coordinate dilation, 
establishes a unitary equivalence between the two Schrodinger operators 

v^/(M+^)(-dydq' + U,iq)) and ^/=' -dVdg^ + V;(g), (3) 
where ^ ^ g-{M+mN+2) y^^^^ '^^ ^ ^^m _ 

Thus, eq.(j2I) is equivalent to 

(^-^ +qN + ^qM ^ ^j^(^) ^ ^ ^ ^_(M+2)/(iV+2)^ ^ _ _^2/(M+2)^^ (5) 

In this transformed Schrodinger equation, the interaction term is now vq^^ and has 
the lower degree, so that v can act as a regular deformation parameter; the former per- 
turbative regime g ^ translates as the asymptotic v —>■ +00 regime. However, at no 
finite v is the problem (jS)) solvable in any traditional sense, and this has severely limited 
practical uses of this reparametrization. On the other hand, this deformation can be fully 
studied by exact WKB analysis, which now handles general (ID) polynomial potentials. 
One earlier detailed study of this sort is based on resurgence theory j2j. Another such 
path from exact WKB to perturbation theory lies in proving the Zinn- Justin conjectures 
about multi-instantons jSiE]- Here, continuing a different type of study initiated in [7j 
(Sees. 3-4) (within an exact WKB framework built upon Sibuya's formalism [H]), we 
seek to specify how the spectral determinants themselves (and related spectral functions) 
asymptotically depend on the coupling parameter v — > +00 (or g ^ 0). Spectral func- 
tions being symmetric functions of all eigenvalues Ek{g) together, the non- uniformity in k 
of perturbative approximations must show up somehow, and the g behavior of such 
objects might not be obviously traceable to existing (fixed-fc) perturbative results. 

As a very concrete example, we may ask: how do the spectral zeta functions Zg{s) = 
J2'kLo^k{g)~'^ precisely behave for g ^ 07 Specially at s = 1 when M = 2: then, that 
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series converges \/g > but term by term it becomes the divergent (odd) harmonic series 
^^(2fc + at g = 0. The latter admits one fundamental regularization by means of a 
sharp summation cutoff K: 

K-l 

J]^^~|(logK + 7+21og2), i^^+oo, (6) 

k=0 

which is (in just a slight variant form) the basic definition of Euler's constant 7. Now, the 
eigenvalues themselves obey Ek{g) ~ 2fc+l for g {k fixed) by perturbation theory, but 
Ek{g) oc [(7(2/^ + l)^]2/(^+2) for k +00 (g fixed) by the asymptotic Bohr-Sommerfeld 
condition, the crossover zone being roughly located by 2Kg + 1 oc (y(^2/(Af-2) Hence 
the series for Zg{l) is another natural, "soft" regularization of the odd harmonic series, 
beginning to act around k oc Kg. Substituting K = Kg into eq.®, we are led to expect 
Zg{l) ~ — jyzi log5' + ^Af- Such an intuitive approach may work for the logarithmic slope, 
but not for identifying the additive constant C^. By contrast, exact WKB analysis can 
yield a precise asymptotic prediction for this (and other) zeta- values: see our final formula 

dsni). 

The outline of the paper is as follows. Sec. 2 recalls essential prerequisites and def- 
initions for the exact WKB approach to be used here. Sec. 3 presents the asymptotic 
problem and its conceptual resolution by exact WKB theory as eq. (I34j) . Sec. 4 performs 
the key computational steps: a class of specific improper action integrals Jq°° n(g) dg 
(where Il{q) are essentially classical momentum functions, and the integrals are primi- 
tively very divergent) are explicitly evaluated. Finally, Sec. 5 processes all intermediate 
calculations into concrete formulae, mainly eqs. (j3^ . (j3Hj) . 



2 Some notions from exact WKB theory 

We recall the essential facts and notations to be used later concerning the exact WKB 
treatment of Schrodinger operators on L^(M), of the form H =^ — d^/dg^ + V{q), with 
a polynomial potential V{q) = + [lower-degree terms], here taken real and even. 
Details and justifications have to be omitted (cf. 7] Sec. 1, and references therein). Such 
operators are self-adjoint, have compact resolvents and commute with the parity operator 
[q I— > —q). A frequently needed quantity (which we call the order of the problem) is 

Miv) = 1 + ^- (7) 

As standard notations, we will also use iplz) = [r' / r](-2), and 7 = Euler's constant. 

2.1 Improper action integral, and residue function 

Important quantities enter at the classical dynamical level around the (complexified) 
momentum function, 

T^xiq) = (Viq) + \f'\ (8) 
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where the constant (—A) stands for the classical energy and, say, A > — inf V (initially). 
Next, improper action integrals over semi- infinite paths prove very useful: J^°° ^xi<l') dg' 
(primitively divergent) is very naturally redefined as the analytical continuation to s = 0, 
when this is Gnite, of 

r+oo 

Ig(s, A) = / {V{q') + A)-'+^/2 dq' (convergent for Re (s) > fi{N)). (9) 

Jq 

Now, the s-plane singularities of Iq{s,X) entirely come from the large-g behavior of the 
integrand. Specifically, the q — * +00 expansion (explicitly computable order by order) 

(Viq) + A)-+^/2 ~ Yl M')<i'"''' (P = ^/2, N/2 - 1, ■ ■ ■) (10) 

p 

implies the singular decomposition 

p+l-Ns 

/,(.,A)^-^/?,(.) ^^_^^ ; (11) 
p 

hence at s = 0, /^(s. A) has at most a simple pole, generated by the p = — 1 term (if any): 

Res,=o/,(s,A) =/?_i(0)/iV, (12) 

a value actually independent of A (save when N = 2) and of q. 

A central distinction sets in at this point: if the "residue function" /3_i(s) = 0, 
the Schrodinger problem {H + A)^ = will behave more simply ("normal" type, N); 
otherwise, "anomaly" corrections will enter (type A). Wholly generic polynomials {V{q) + 
A) are of type A; still, in a sense, "a majority" of them have f3-i{s) = 0: among the even 
ones, already all those having a degree N multiple of 4 (and, among the non-even ones, 
all those of odd degree). 

Thus, for the N type, f^°° ^xil') dg' can be readily defined as the analytical contin- 
uation of Ig{s, A) to the (regular) point s = 0; whereas for the general (A) type, the best 
specification is not the bare finite part of Ig{s, A) at the pole s = (denoted FP5=o-^(}('5, A)), 
but rather ([7 , eq.(32)) 

r+00 

/ Ux{q') dq' = FP,=o/,(s, A) + 2(1 - log 2) /?_i(0)/iV (13) 

Jq 

in order to preserve the basic identities ()21|) below. Additivity is also maintained: 

p+oo pq" p+oo 

/ Ylxiq') dq' = / Uxiq') dq' + / Uxiq') dq' for all g, q" finite (14) 

Jq Jq Jq" 

(because a finite integral J^^ (^(o'O + A)^^^^/^dg' is entire in s). 

Remarks: (i) J^°° ^xW) dq' is an "Agmon distance from q to +cxd", suitably renormal- 
ized; (ii) this procedure is a classical counterpart to zeta-regularization at the quantum 
level; (iii) just as the extra term in eq. ()13p . all anomaly terms here will simply be propor- 
tional to the value /9_i(0), but more general forms occur elsewhere [7]. 
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2.2 Spectral functions 

An operator H as above has a purely discrete real spectrum {Aq < Ai < A2 < ■ ■ ■}, 
(Afc t +00), where even (resp. odd) k correspond to eigenfunctions of even (resp. odd) 
parity. Parity-split spectral zeta functions (a la Hurwitz) can be defined as 

A) = (Afc + A)-^ for Re s > fi{N) (15) 

I, even 
odd 

and, say, A > — Aq; many results however take a sharper form upon a skew versus a full 
zeta function, respectively defined as 

Z^ = Z+-Z-, Z = Z+ + Z-. (16) 

Spectral determinants -D^(A) = det^{H + A) are defined by zeta-regularization, as 

D^{\) = exp[-d,Z^is,\)],=o {and = D+ / D- , D = D+D-), (17) 

where s = is reached by analytical continuation from {Re s > fi{N)}. More constructive 
specifications are 

- their Weierstrass infinite products (written for fi{N) < 2, which is true here): 

D^{\) = D±(0)e^P»=i^^('^'°)^ Y[ (l + A/Afc)e-^/^^ (18) 

■ even 
odd 

= D^{0) Ylil + A/Afc) when /i(iV) < 1, i.e., > 2 (19) 

I even 
odd 

and likewise for D, D^; this shows that the determinants continue to entire functions (of 
order fi{N)) in the variable A (except for D^, meromorphic); 

- the basic identities of the exact- WKB method: let \1'a(q') be the canonical recessive 
solution of the differential equation, specified through its q +00 asymptotic form 

^A(g)-n,(g)-V2e/rn.(.')d.'^ (20) 
where J^"^ dq' is fixed according to eq. lfl^) : then, under that precise normalization, 

D-{X) ^ ^,{0), D+{X) ^ -vl/',(0), (21) 

(also valid for a rescaled potential, i.e., V{q) = uq^ + ■ • •, with u > 0). Remark: the 
solutions obeying are close to Sibuya's subdominant solutions but the two nor- 
malizations fully coincide only when the type is N. 

Finally, we will need the transformation rules for spectral functions under a global 
spectral dilation (A^ t— > rA^, r = est. > 0). Obviously, 

Z^{s, A) ^ r-'Z^{s, A/r) for Re s > ij{N) (22) 

(and likewise for Z, Z^); hence upon continuation to s = 0, and applying eq. p7|) . 

D{X) ^ r^(°'^/'-)D(A/r), ^^(A) ^ r^^^^'^/'^D''{X/r) (23) 

where, moreover, ([7], eqs.(27),(37)) 

Z{0, A) = -2(3_,{0)/N, Z^{0, A) = 1/2. (24) 
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3 The asymptotic v +(X) problem 

We now return to the Schrodinger operator H{v) = — d^/dg^ + + vq^\ as in eq.© 
(A^ > M > 2 both even, f > 0). We will find the asymptotic behaviors of its spectral 
determinants in the regime of singular perturbation theory for eq.(|H): 

D^{\ v) = det^{H{v) + A) (A > 0), v = ^ (25) 

To lowest-order in the g ^ perturbation theory, the individual eigenvalues Afc(f) of 
H{v) become asymptotic to those of Ho{v) = —d^/dq^+vq^. We then expect <let^{H{y) + 
A) to become somehow asymptotically proportional to dei^ {Hoiy) + A) as f — +oo, but 
the latter regime is singular and moreover non- uniform in k\ hence the actual behavior of 
the determinants cannot be taken for granted. In 7] (Sees. 3-4), we tackled it for a few 
binomial potentials and exclusively at A = 0; now we will do it in full generality. 



3.1 Detailed anomaly types 

As argued in Sec. 2.1, it is essential to distinguish between normal (zero-residue) and 
anomalous (non-zero residue) cases, but this now applies independently to the coupled 
(= H{y)) and the uncoupled (= Hoiy)) problems. 

- the coupled problem (Jlx{qY = q^ + vq^ + A): the residue function is the coefficient 
of g-i-^* in the generalized binomial expansion for g^(i/2-'5)(i -|- vq^^^'^ + Ag^^)^/^~*. 
When > 2 as here, the residue function cannot involve A; specifically, 

P-i{s) = unless ~*" = j eW ("anomaly condition of level j"); (26) 

thus, anomalies attach to specially correlated exponents N, M only: 

(level j :) A^ = 2jm — 2, M = N — m for m G N* (with m even for even potentials); 

(27) 

and then 

- the uncoupled problem (Jl^ x^qY = vq^ + A): the same calculation now simply yields 

(3_^{s) = v-^/^X{l/2 - s) ifAf = 2 [AiforA^O], otherwise /3_i(s) = [N]. 

(29) 

The harmonic oscillator (n(g)^ = vq"^ + A) thus gives the prime example of anomaly, 
actually the unique case (among all potentials) where the residue depends on the spectral 
parameter; all other binomials {vq'^^ + A} (M ^ 2) are of type N. 

The type can abruptly change either way in the v +oo limit, giving birth to four 
distinct variants (the "basic" example of eq.(IT} is not the simplest!): 
N ^ N: e.g., V{q) = q^ + vq^; 
Aj N: e.g., V{q) = q^ + vq^, of level j = 2; 
N Ai: e.g., V{q) = q'^ + vq^ (the "basic" example) when A 7^ 0; 
Aj —>■ Ai'. only one case, V{q) = q^ + vq^ when A 7^ 0, for which j = 1. 
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3.2 The main estimate 

We can relate the coupled and uncoupled spectral determinants very easily through a key 
result of exact WKB theory, the basic identities (j2U)- These are to be written for both 
(coupled and uncoupled) problems independently: 

det-iHiv) + \) ^ ^x{0,v), det+(if(t;) + A) = -^^0,^^), .^^1 

det-(^o(t^) + A) = ^o,A(0,t;), det-^{Ho{v) + X) = -%^^{0,v), ^ ' 

where '^x{q,v), resp. \E'o,a(5', "w) are the canonical recessive solutions of {H{y) + A)\E' = 0, 
resp. {Ho{v) + A)\E'o = 0. So, the problem boils down to relating '^x{q,v) and \E'o,a(?, "w) 
near g = as f ^ +oo. 

Now, as soon as ^ v, the term becomes a negligible perturbation of vq^ 

within the Schrodinger equation, hence the recessive solution \E'A(q', f) has to become 
asymptotically proportional to \E'o,a('?, "w) in that regime (given that the WKB form ()2n|l 
holds asymptotically for Il\{q) +oo whichever way the limit takes place, including 
V +00 at fixed q). The only problem is then to determine the asymptotic ratio 
'^x{q,v)/'^o^x{q,v) ~ C{X,v) (g-independent) as v ^ +oo at fixed q. By contrast, the 
alternative normalization of recessive solutions based at q = 0, 

^x{q,v) ^ Ux{q,v)-^/^ e- ^o^'i'i''-)^i\ % x{q,v) Uq x{q,v)-^/'' e-^o^'Mi'^-)dg' 

(31) 

(for q +oo) immediately entails 

^xiq, v) ~ %,xiq, v) for v +oo, \qf-'' < v, (32) 

because the asymptotic equivalence Ilx{q', v) ~ no,A(?', v) can be used all over the bounded 
interval [0, q]. 

The final issue is to relate the two normalizations, the canonical one of eq. ()2U|) ( "based 
at g = +oo") and the latter one based at g = 0. Thanks to eq. lfT^ . the answer is simply 

^x{q, v) = ei'o'" n,(g',^)dg' ^^^^^ ^^^^ likewise for ^o,a with Uo,x). (33) 

Finally, putting together eqs. ()30|) - (jn^ . we end up with the comparison formula 

det^(^(t;)+A)~e^(^'^)e-^o(^'^)det±(^o(^^) + A) {v ^ +oo) (34) 

(stated in most general terms), where 



SiX,v) = / nA(g,f)dg, resp. So{X,v) = / no,A(g,f)dg, (35) 
Jo Jo 

are coupled, resp. uncoupled, improper action integrals. Specifically here. 



SiX,v) 



/• + 00 p + oo 

/ (g^ + vq^' + A)^/2 ^ggp^ g^^x, v) = / {vq^' + A)^/^ (gg) 

Jo Jo 



The problem has thus been decomposed and reduced to the asymptotic {v +oo) eval- 
uation of the two action integrals of eq. ljHBj) separately. 
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4 Explicit formulae for improper action integrals 



This Section constitutes a kind of technical digression, but the effective computations 
of improper action integrals to be presented might also be of autonomous interest and 
applicability. 



4.1 Binomial 11 (g)^ : exact evaluation 

We compute the improper action integral n(g) dg exactly for a binomial Il{q)^ = 
uq^ + f g*^, in the rather general setting A^>M>0, m, f>0, resulting in the multi- 
purpose formulae (pn|) and (jl^ (here, and M might even be non-integers). 
At the core, j^'^ n(g) dg = lims=o Io{s) where 

/•+00 

Io{s) = / (uq'^ + vq'Y'-'dq {Re s > ^ + i), (37) 
Jo 

as long as the limit (understood as the analytical continuation to s = 0) is finite. Now 
the right-hand side reduces to a Eulerian integral, of the form 

r+oo 

/ (ax + 6)^/2""a;°-Ma; = a-"6i/2+"-"r(a)r(s-a-l/2)/r(s-l/2) (38) 
Jo 

(under the change of variable q^''^^ = u~^vx] here, a = [M(l — 2s) + 2]/[2(A^ — M)]); 
more precisely, 

„/M(l-2£)+2x Af(l-2^)+2 ^ 

^ ' {N -M) r(s-l/2) ^ ' 

Consequently, at s = 0, 

/.+00 p/ Af+2 ^1 p/ N+2 \ 

Jo (Ar-M)r(-l/2) ^ ' 

in the normal case, i.e., when the right-hand side is finite, meaning here 2(jV"~^M) ^ ^' 
As concrete examples of this N type: 

r+oo 

/ {q' + vq'y/^dq = -v'/y3 (41) 
Jo 

r+oo 

/ (ug^ + A)^/Mg = -(2v/^)-Y(l + i)r(-|-^)«-^AH^ (AT ^2). (42) 

Now, the right-hand side of eq. ()40|) turns infinite whenever (2j — 1)A^ = 2{jM + 1) 
for some j G N* (j = cannot occur); this is precisely the anomaly condition Aj of level 
j. The binomials of any type Aj can be readily (albeit tediously) handled by applying 
eq. lfT^ to Io{s). (The cases with j = 1 as well as the action integral of eq. pTj) were 
implicitly evaluated in jZl, by a different route.) First, the residue is 
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then, the finite part at s = of eq. (jH^ gets extracted as 



/5-i(0) 



2j 



N+2 



M / 

^(j + 1) - \ogv + — (^-^(j - 1/2) + \ogu 



N 



^(-1/2) 



so that finally, when M = [{j — 1/2)N — l]/j (condition A^), eq. (fT!?|l yields 



[uq +vq ) ' dq= 



1 2M / , 4^ 1 

logt; + ^- + — (log2+2logM-^ 



m=l 



m=l 



2m- 1 



(44) 

The cases with j = 1 are of special interest. Besides the harmonic oscillator, the 
general binomials of type Ai are just the supersymmetric potentials (at zero energy): 



U{qY = uq^ + vq^ with N = 2M + 2 (M > 0), 
and eq.(jlH) distinctly simplifies to 



(45) 



+ 00 



{uq^+vq^/'-y/^dq 



N — 2 

_log^; + l + __(log2 + ilogn) (j = l). (46) 



In particular, for the harmonic oscillator (A^ = 2) at a general energy value (—A), 

DO 

(V + xy/^ dq = ii;-'/'A(l - log A). (47) 



4.2 Trinomial Il{qy : asymptotic v ^ oo evaluation 

We now consider a trinomial n(g)^ of the form q^ + vq^^ + A, with even > M > 0, and 
a systematically constant third term: the spectral parameter itself, A (> 0) (= minus the 
total energy). One of the coefficients can always be scaled out to unity, and we have done 
this for the highest power initially. 

In this case we can no longer compute the action integral J^°° 11 (g) dg exactly. In 
view of eq. ()34|) . however, we mainly need its large-f behavior, specially for v +oo in 
order to recover singular perturbation theory according to eq.Q (but as in [7j, we expect 
the results to remain valid over suitable sectors in the complex t;-plane). 

According to the zeta-regularization idea, we must start from the large-f behavior of 

J(s; A, v) =^ {q^ + vq^ + dg; this problem is rather delicate, so any brute-force 

expansion scheme is dubious. Instead, we apply the following general idea: if the function 
I{v) under study is an inverse Mellin transform. 



I{v) = (27ri)-i / i{a)v''da, 

J c— ioo 



(48) 



then the singularities of I{a) in the half-plane {Re a < c} encode the large-w behav- 
ior of I{v). Thus (by the residue calculus) any polar part of the form A{a — o"o)~^ + 
-B(cr — ctq)^^ in represents an asymptotic contribution v'^°{A\ogv + B) to /(f). This 
is particularly valuable for /(s;A,f), because its direct Mellin transform /(s; A, a) =^ 
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AIM 



-IIM 



-3IM 



■ confluence if M=2 



1/2 
s 



Figure 1: Schematic plot (using some non-integer A^, M) of the poles cr(s) of the Mellin trans- 
form I{s;\,a) in ea.(|49|). The two main contributing poles (in the v — > +00 limit) are drawn 
with bold lines; non-contributing poles are drawn with dashed lines. 



J(s; a, v) v'"''^ dv is exactly computable (using the same formula as for the ex- 
act action integral of a binomial but now twice in succession), and it is a meromorphic 
function of cr: formally, 



dvv 



(g^^+V^+A) 



a/2- 



r -a r s + a- 



r(s-i/2) 



I(s: X,a) 



r -fx r 



■ Mcr+1 - 

, AT , 



r(s+(T-|- 



Mcr+1 ' 

N ■ 



A" 



(49) 



A^r(s-l/2) 

(using the change of variable = Xr for the g- integration). Now this Mellin transform 
also has to genuinely exist somewhere; here, all integrations converge in some strip a' < 
Re 0" < provided Re s > n{N), and the inverse transformation ()48|) applies with c = — 0. 
Consequently, the poles a{s) relevant to the current asymptotic problem are those which 
lie in {Re a < 0} when Re s > n{N), and their contributions are then to be analytically 
continued to s = 0. Overall, the poles in eg . (j49|) form three arithmetic progressions, 
one for each Gamma factor in numerator; they are real for real s (fig. 1). At s = 0, 
any pole a{s) will contribute an asymptotic term of degree = a{0) in v (on general 
grounds) and of degree dx = ^ + — ^^^^ cr(0) in A (by examination of eq. pUj) ). At the 
end, we plan to keep the terms of degree dg < in the perturbative coupling constant g 
(discarding o(l) terms when g 0); now the Symanzik scaling (eq.(jSI) at fixed E) entails 
dg = —^^{dv + jf^dx) = — (Mcr(O) + 1)/A^; altogether, dg <0 then amounts to keeping 
only the poles for which a{0) > —1/M. 

When M > 2 (as here), only two poles a{s) satisfy both criteria, (in real form) a{s) < 
for s > /i(iV) and a{0) > —1/M: they are, in decreasing order (at s = 0), 



N 



ao[s) = 



N-M 



1 1 

2 + iV-^ 



(leading) and Cils) = — — (subleading). (50) 



M 



They are generically simple, with two exceptions at s = 0: o"o(0) = 2{N-l4) becomes 
confluent with the (fixed) pole a = +j when the coupled problem is of anomalous type 
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Aj] and independently, ai becomes confluent with the (next mobile) pole j^^j^ji^ ~^ J/ ~ 
s — 1) when the uncoupled problem is anomalous, i.e., M = 2. The latter will induce a 
usual double-pole contribution; the former confluence is worse, making the inverse-Mellin 
representation singular as the integration path gets pinched between the two confluent 
poles. 

We now specifically evaluate the two dominant polar contributions, from cxo and cxi. 

- the leading pole cro(O) = 2(n-'m) • coupled problem is of type N this pole remains 
simple, and its asymptotic contribution [ReSo-Q /(s; A, cr)] t>°"° turns out (by inspection) to 
be just Jf^'^iq'^ + vq'^^Y^'^ dg (as given by eq. pUj) aX u = 1). Furthermore, d\I{s] A, a) oc 
I{s + 1; A, 0"), an operation which precisely annihilates this leading pole part in all cases, 
so the latter has to be a constant in A; then, its computation at A = precisely yields 
/g^°°(g^ + vq^'^Y^"^ dg, now including the conHuent cases (Aj). 

- the subleading pole ai = —1/M: if M > 2 this pole remains simple, and its asymptotic 
contribution [ReSg-i /(s; A, cr)] w'^^ coincides with J(^°°{vq^^ + A)^/^ dg as given by eq. ip^ . 
Under confluence (M = 2), the contribution becomes that of the double pole of eq.(|l^ 
at (T = —1/2: this is jQ°°{vq^ + A)^/^ dg + -j^^Ai{\, v) where the action integral is given 
by eq. ljTTj) . and 

A^{\v) = it;~^/2^(logt; + 21og2). (51) 
All in all, the asymptotic v ^ +oo formula for the trinomial action integral is then 

f+OO f + OO f + OO AT 

/ (g^ W+A)^/Mg ~ / (g^W)'/'dg+/ (vq^'+Xy/' dq+SM,2 ^ A^iX,v), 
Jo Jo Jo i\—z 

(52) 

where °°(g^ + vq^^Y^"^ dg is specified through eq. pUj) if the coupled problem is of type 
N, or else eq. p^ if the coupled problem is of type A^ (i.e., if 2(jv-M) ~ J ^ ^*)! whereas 

Io°^in^' + Xy/^dq is given by eq.^ if M > 2, or eq.(|lZj) if M = 2 (and Sm,2 is a 
Kronecker delta symbol). 

5 Asymptotic behaviors of spectral functions 

The theoretical results of Sees. 3-4 translate into concrete formulae for spectral functions 
in the v +oo regime. 

5.1 The spectral determinants 

Upon substituting the explicit formulae of Sec. 4 into eq.(j21I), /g^°°(fg^ + A)^/^ dg cancels 
out, and a slightly simpler v —>■ +oo formula results: 

det±(-dVdg^+g^+ V+A) ~ e^^°°(^"+^''")^'''^^+^-^^ iv^^i(a,^) det±(-dVdg2+ V'+A), 

(53) 

where Jg °°(g^ + vq^Y^^ dg is given through eq. ljlUj) if the coupled problem is of type N, 
or (jUj) if it is of type A^, and Ai{X,v) by eq.(jSH). 

Being homogeneous, the uncoupled potentials obey a simpler form of the scaling eq.©: 
(— d^/dg^ + fg*^) is unitarily equivalent to v'^f'^^^'^\—d^/dq^ + g''^); then the scaling laws 
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()23|) apply with r = ti2/(^'-f+2)- these laws are more awkward for the det''^ than for the 
full and skew determinants det and det^, we now switch to the latter combinations and 
explicitly get 

det(-dVdg' + vq"" + X) = det(-dVdg2 + g*^ + t;-2/(A^+2)A) (M ^ 2) (54) 

det(-dVdg2 + vq^ + X) = t;--"'^'V4 det(-dVdg2 + q^ + v'^/^X) (55) 
detP(-dVdg2 + vq^' + X) = v'/'^^'+''> de^i-dydq' + q'' + v-^'^'^'^^^ X). (56) 

Remark: the harmonic-oscillator determinants are actually known in closed form ([2j, 
eqs.(155)); e.g., eq. (|HH|l has the fully explicit form (to be used later) 

det(-dVdg' + V + A) = i;-^"''"^/^2-^"'^'^/2^/ p( i ^ v'^'^X)). (57) 

Thus, eqs.dSS} plus supply the v — » +oo behaviors at fixed A of the coupled 

determinants in terms of the corresponding uncoupled determinants at A = (which are 
computable numbers, cf. [Zj, eq.(136)). 

However, our main concern is rather the singular perturbation limit: v = g~(^^+'^y(^+'^^ 
+00 with t)~2/{M+2)j^ 4|f (^_^) fixed, according to eq.©. The exphcit final results, de- 
duced from eqs. lj^Hjl - lj^Uj) after rescaling both sides, are then 

det(-dVdg' + q^' + gq"" - E) / det(-dVdg' + q^' - E) 

^ ^-|wTle2/o^°°('?''+''''")'^''i'?e-^«-2i^(^'°s''+^'°*52)s for type A,- 

(/3_i(0) (given by eq. ljlHj) ). and type, both refer to the coupled problem); whereas the 
skew determinants always behave straightforwardly: 

detP(-dVdg2 + g*^ + gq^ - E) ^ detP(-dVdg2 + g^'' - E). (59) 

The main result here is the explicit non-trivial prefactor in eq. (j58p . Its essential singularity 
for g should relate to the non-uniformity of this limit with respect to the quantum 
number k. By contrast, the dependence of its logarithm upon E is elementary, consisting 
only of (a) constant term(s) (already determined in 17] for some cases), then a linear term, 
and nothing else. The basic example (^, being of type N, thus gives 

det(-dVdg' + + gq^ - E) ^ e'^/^^ e(i°g9/2-2iog2)£; ^g^^_^2/^^2 ^ ^2 _ (gg^ 
(Note the "instanton-like" structure of the first prefactor, computed by eq. pTjl .) 

5.2 The spectral zeta functions 

Over the spectrum {Ek{g)} of the rescaled operator {—d'^/dq'^+q'^^+gq'^), we can consider 
the full and skew spectral zeta functions 

00 00 
Z,{s- E) J2^E,ig) - E)-^, Zj(s; E) 5^(-l)^(i?,(^?) - E)-^ (61) 

fc=0 fc=0 
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Figure 2: Illustration of the two g ^ behaviors of spectral zeta functions at E = 0, Zg{s) = 
Zg{s; 0) and Z^{s) '= Z^{s; 0), as computed from numerical spectra (by means of eas. (|HT|) ') for a 
sample of potentials V{q) = +gq^ with A'^ = 4 (+) or 6 (x). Left: regular behaviors verifying 
eq. (|63j) . plotted against \/v (the variable which empirically appears to yield the simplest depen- 
dence); remark: Z^(l) = 7r/4, Zo(2) = vrVS, and Z^(2) = Catalan's constant (ss 0.9159656). 
Right: singular behavior of Zg{l), plotted in a log 5 coordinate; the straight lines show the 
theoretical asymptotic predictions (ea. H66() l. 



for, say, integer s G N*, in which case they converge for g > and relate to the spectral 
determinants in a simpler way than for general s, 

Zgis; E) ^ -^-^^ log det(-dVdg^ + q'' + gq'^ - E), (62) 

(obtained from eq. lfT^ upon rescaling; and likewise for (Z^, det^)). 

Assuming all previous estimates are stable under ii^-differentiations (as is usually the 
case in WKB theory), the preceding formulae imply the regular behaviors (see fig. 2, left) 

Zg{s-E)r~.Z,{s-E), Z^{s;E)r^Z^{s;E) (g ^ 0), (63) 

except for Zg{l;E) (the resolvent trace) when M = 2, which gives the singular case 
{Zq{1;E) infinite, while Zq{1;E) stays finite). Those patterns were conjectured in [7j 
(Sec. 3), but not the precise divergent behavior of Zg{l;E), which requires to know the 
ii^-linear term in the exponent of the determinant ratio ()58|1 . For s = 1, eq. ljH^ needs to 
be regularized at (7 = 0, as 

- (d/dE) logdet(-dVdg' + q^ - E) = -|[^(|(1 - E)) + log2] (64) 

(using the known closed form (jSZj) of the harmonic-oscillator determinant). Then the 
logarithmic differentiation of eqs. (j58|) for M = 2 yields the g ^ behavior of Zg{l;E) 
for all potentials + gq^ (irrespective of type), as the following singular expression: 

Zg{l; E) ~ \ogg + N\og2) - \ [^(i(l - E)) + log 2]. (65) 
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For instance, at £" = this gives (see fig. 2, right) 

J^Ekigr' ~ -^^log^? + -(7+^^^log2) (g^O) (66) 

fc=0 

ilog5(+ i(7+51og2) forA^ = 4, -i loggf + ^7 +21og2 for iV = 6, . . . 

(to be compared with the sharp cutoff regularization of eq.®)- 

5.3 Concluding remarks 

We have completed here one "exercise in exact quantization" begun in |7|: we gave the 
g ^ behavior of the spectral determinants det^(— d^/dg^ + q^^ + gq^ — E), now for 
general parameter values. While it may appear wasteful to use a wholly exact approach 
for perturbative calculations, exact WKB analysis actually proved quite efficient for the 
task; inversely, such problems help to strengthen the practical sides of that field, which 
still need further development. 

We are also confident that the above approach can be extended further, both to 
complex parameter asymptotics and towards higher orders in powers of g. 
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